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Abstract: Many theories, from Supersymmetry to models of Strong Electroweak Sym- 
metry Breaking, look at the production of four top quarks as an interesting channel to 
evidentiate signals of new physics beyond the Standard Model. The production of four-top 
final states requires large partonic energies, above the 4mj threshold, that are available 
at the CERN Large Hadron Collider and will become more and more accessible with in- 
creasing energy and luminosity of the proton beams. A good theoretical control on the 
Standard Model background is a fundamental prerequisite for a correct interpretation of 
the possible signals of new physics that may arise in this channel. In this paper we report 
on the calculation of the next-to-leading order QCD corrections to the Standard Model 
process pp —> ttti + X. As it is customary for such studies, we present results for both 
integrated and differential cross sections. A judicious choice of a dynamical scale allows us 
to obtain nearly constant /C-factors in most distributions. 
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1 Introduction 

Since its discovery in 1995 by the CDF and DO experiments at Fermilab, the top quark, 
the most massive of all the observed elementary particles, has been extensively scrutinized. 
Several properties and observables have been already analysed including, among others, 
the top-pair and single-top production cross sections together with the measurement of the 
top quark mass. Most of these measurements have been, however, limited by the relatively 
small statistics collected at the Tevatron. With a much higher center of mass energy and 
luminosity, the Large Hadron Collider (LHC) has all the features of a top factory and 
will significantly add to the previous measurements in this field. At the nominal design 
value of the beam energy (-y/s = 14 TeV) and assuming a luminosity of 10 fb _1 /year, the 
LHC will collect about 8 million top anti-top pair events and another few million tops 
via electroweak single top quark production. In fact even at present, running at ^/s = 7 
TeV, it has made possible to start re-examining the main properties of the top quark with 
significant precision. 

The large statistics available will open a window on entirely new measurements as well 
as on the analysis of novel, more complex final states. As to titt production, the LHC 
energy is sufficient to produce such events at a sensible rate. Assuming only the standard, 
dominant t — > Wb decay, this process leads to W + W~W + W~bbbb final states with leptonic 
and/or hadronic decays of gauge bosons. The experimental signature of such final states 
is characterized by the presence of opposite and same-sign charged leptons, large missing 
energy and jets with a considerable number of b-jets. These signatures are precisely the ones 
that are expected in the Higgs boson(s) and new physics searches. Thus, titi production 
is an interesting channel to probe several realizations of Standard Model extensions at the 
LHC. Some prominent examples include models of Higgs and top compositeness as well 
as models involving the production of new colored resonances with large couplings to the 
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top quark. In the latter case, the most studied ones are Kaluza-Klein gluons from the 
Randall-Sudrum warped extra dimensions. Many of these models predict effective four-top 
quark interactions as well as new processes such as pp — > GG and pp — > tt G, where G is a 
new heavy particle decaying to ti, leading to titi final states (see e.g. [1-10] and references 
therein). In this context, a precise theoretical description of the four-top production rate 
in the Standard Model may help to constrain new physics scenarios. 

In addition, titi is a major background for many processes arising from supersymmet- 
ric extensions of the Standard Model, the most noticeable example being the production 
of a heavy Higgs boson [11, 12]. An accurate theoretical description of four-top production 
would help in this case to determine the Higgs boson self coupling and thus contribute to 
better understand the Higgs boson potential in this minimal SM extension. Furthermore, 
multi-top final states can also be produced via long cascade decays of colored superymmet- 
ric particles, such as squarks or gluinos [13, 14]. 

Last but not least, titi production is the last process in the so-called Les Houches next- 
to-leading order (NLO) experimenter's wishlist [15] that is still waiting for a calculation. 
Recent breakthrough in one-loop calculational techniques, sometimes referred to as "the 
unitary revolution" [16-19], together with great improvement in more traditional methods, 
have led to a tremendous progress in the calculation of multi-leg processes at hadron 
colliders. This is exemplified by the calculation of the following 2 — > 4 processes: pp{pp) — > 
tibb+X [20-23], pp(pp) -> tijj+X [24, 25], pp(pp) -> W + W~bb+X [26, 27], pp -> bbbb+X 
[28], pp(pp) -> W+W-jj + X [29, 30], pp -»• W + W + jj + X [31, 32], pp(pp) -»■ W + 3j + X 
[33, 34], pp -»• Z/ 7 * + 3j + X [35], pp -> W'yyj + X [36] and pp -»• 4 j + X [37]. In 
addition, the first NLO QCD corrections to 2 — > 5 processes, i.e. pp — > W + 4j + X and 
pp — > Z/^f* + Aj + X have recently been completed [38, 39]. 

In this paper we present the results of a calculation of pp — > titi + X in the Standard 
Model at NLO QCD accuracy. This calculation makes the Les Houches wishlist finally 
complete, with at least one independent calculation available for each benchmark process. 
The paper is organized as follows. In Section 2 we briefly describe the details of our 
calculation. Numerical results for the integrated and differential cross sections are presented 
in Section 3. Finally, in Section 4 we give our conclusions. 

2 Details of the calculation 

At the leading order (LO) in perturbative expansion, titi final states are produced via the 
scattering of either two gluons or one quark and the corresponding anti-quark. A represen- 
tative set of Feynman diagrams contributing to the pp — > titi process at 0(aj.) is depicted 
in Figure 1. In total, there are 72 LO diagrams for gg — > titi and 14 for qq — > titi. Even 
though we do not actually employ Feynman diagrams, it is customary to present them as 
a measure of the complexity of the calculation. The calculation of scattering amplitudes is 
based on well-known off-shell iterative algorithms [40-42] , performed automatically within 
the Helac-Dipoles package [43] and cross checked with the Helac-Phegas Monte Carlo 
program [44-46]. A perfect agreement has been found in all cases. Phase-space optimiza- 
tion and integration has been performed with the help of Parni [47] and Kaleu [48]. 
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Figure 1. A representative set of Feynman diagrams contributing to the pp — > tttt process at 0{a A s ). 
Double lines correspond to top quarks, single lines to light quarks and wiggly ones to gluons. 




Figure 2. A representative set of pentagon and hexagon diagrams for the pp — > tttt + X process 
at NLO QCD. Double lines correspond to top quarks, single lines to light quarks and wiggly ones to 
gluons. 




Figure 3. A representative set of Feynman diagrams contributing to the real emission corrections 
to the pp — > tttt + X process at O(a^). Double lines correspond to top quarks, single lines to light 
quarks and wiggly ones to gluons. 



At the NLO level, the virtual corrections are obtained from the interference of the one- 
loop diagrams with the tree level amplitude. They can be classified into self-energy, vertex, 
box-type, pentagon-type and hexagon-type corrections. A representative set of pentagon 
and hexagon diagrams is shown in Figure 2. The number of one-loop Feynman diagrams 
for the pp — > tttt process, as obtained with Qgraph [49], is 2200 for gg — > tttt and 410 for 
the qq — > tttt. Virtual corrections are evaluated in d = 4 — 2e dimensions in the 't Hooft- 
Veltman version of the dimensional regularization scheme. The singularities coming from 
infrared divergent pieces are canceled by the corresponding ones arising from the so-called 
integrated dipoles, i.e. from the counterterms of the adopted subtraction scheme integrated 
over the phase space of the unresolved parton. The finite contributions of the loop diagrams 
are evaluated numerically in d = 4 dimension. Finally, to ensure numerical stability, we 
check Ward identities for every phase space point. The calculation of the virtual corrections 
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Partontc; 
SUBPROCESS 


TMttmrfr Of Feynman 
Diagrams 


Nttmrfr Of 

J. i| KJ 1\L J_J J_J J. V V/ 1 

DlPOLES 


gg -)• tttig 


682 


30 


qq — > ttttg 


128 


30 


gq — > ttttq 


128 


10 


qg — > ttttq 


128 


10 



Table 1. Partonic subprocesses contributing to the subtracted real emissions at O(a^) for the 
pp — > titi + X process. Also, the number of Feynman diagrams and the number of Catani- Seymour 
dipoles corresponding to these subprocesses are presented. 

is achieved with the help of the package Helac-ILoop [50] which incorporates CutTools 
[51, 52] and OneLOop [55] as cornerstones. The first code contains an implementation of 
the OPP method for the reduction of one-loop amplitudes at the integrand level, while the 
second one is dedicated to the evaluation of the one- loop scalar functions. Renormalization 
is done, as usual, by evaluating tree-level diagrams with counterterms. For our process, we 
choose to renormalize the coupling in the MS scheme with five active flavors and the top 
quark decoupled. The mass renormalization is performed in the on-shell scheme. 

The real emission corrections to the LO process arise from tree-level amplitudes with 
one additional parton, i.e. an additional gluon, or a quark anti-quark pair replacing a gluon. 
All possible contributions can be classified into the four categories presented in Table 1, 
together with the number of Feynman diagrams and the Catani-Seymour dipoles corre- 
sponding to each subprocess. Typical examples of the real emission graphs are displayed 
in Figure 3. 

For the calculation of the real emission contributions, the package Helac-Dipoles 
[43] is employed. It implements the massless dipole formalism of Catani and Seymour [56], 
as well as its massive version as developed by Catani, Dittmaier, Seymour and Trocsanyi 
[57], for arbitrary helicity eigenstates of the external partons. A phase space restriction on 
the contribution of the dipoles as introduced by Nagy and Trocsanyi [58] is also included. 
For the real corrections as well, we adopt Kaleu equipped with additional, special dipoles 
channels that proved to be important for phase-space optimization in the subtracted real 
emission part [26]. 

To summarize, our computational system relies on Helac-ILoop and Helac-Dipoles 
that are both parts of the Helac-NLO framework for NLO QCD calculations [59]. The 
framework is publicly available 1 . All numerical results are obtained using the same meth- 
ods as presented in our previous work, therefore we will not describe them here in more 
detail. Rather, we address the interested reader to our earlier publications in this field [22- 
26]. We would like to emphasize, however, that all parts of the calculation are performed 
in a completely automatic and fully numerical way. 



http : / /helac-phegas . web . cern . ch/helac-phegas/ 
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Process 


VLO [fb] 


<Eo =1 ^] 


<Eo =aU1 [fb] 


K-Factor 


[%] 


pp -> titi+X 


12.056(6) 


15.33(2) 


15.35(3) 


1.27 


27 



Table 2. Integrated cross section at LO and NLO for pp — > titi + X production at the LHC with 
y/s = 14 TeV. Results for the MSTW2008 PDF set are presented. In the last two columns the K 
factor, defined as the ratio of the NLO cross section to the respective LO result, and NLO corrections 
in % are given. The scale choice is fiF = t l R = A'o = 2mt< 

3 Predictions for the LHC 

We consider the process pp — > titi + X at the LHC with a nominal design center of mass 
energy of \fs = 14 TeV. For the top quark mass, we use the Tevatron average value mt 
= 173.2 GeV as measured by the CDF and DO experiments. Our calculation, like any 
fixed-order one, embodies a residual dependence on the renormalization scale ((Mr) and 
the factorization scale {^f) arising from truncation of the perturbative expansion. As 
a consequence, the value of observables depends on the values of fj,R and fip that are 
provided as input parameters. For many processes, a natural scale can be easily identified 
in several ways. For example, one may consider the mass of the heavy particle appearing 
in the process, or even the typical momentum transfer or the total transverse energy of 
the process. In the case of titi, setting [xq = fiR = fip = 2mt seems a quite natural scale 
choice. In fact, when considering total cross sections, effects of the phase space regions 
close to the top production threshold are expected to dominate, which justifies our choice. 

We consider the MSTW2008 set of parton distribution functions [60] as our default 
PDF set. In particular, we take MSTW20081o68cl PDFs with 1-loop running a s at LO and 
MSTW2008nlo68cl PDFs with 2- loop running a s at NLO, including five active flavors. The 
strong coupling constant is provided by the PDF set itself. It should be noticed that the 
contribution induced by bottom-quark densities amounts to 0.05% at the LO. We decided 
therefore to neglect the contribution of the bb initial state in the computation of the total 
cross section at both LO and NLO. Outgoing top and anti-top quarks are treated on-shell 
without any cut restriction on them. At NLO, an additional final-state parton arises. 
Since this parton has no other partons in the final state to be recombined with (the tops 
are always assumed to be tagged), it will be the only responsible for the possible additional 
jet. As a consequence, the properties of the resulting jet are independent of the choice of 
the jet algorithm and the size of the cone radius. We set no restriction on the kinematics 
of the extra jet. 

3.1 Integrated cross section and its scale dependence for fiQ = 2mt 

We start our presentation of the results with the integrated cross sections. In Table 2 cross 
sections at LO and NLO for pp —> titi + X production at the LHC with yfs = 14 TeV, are 
presented. Our results for /ig = ^ m t scale using the MSTW2008 PDF set are 
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Figure 4. Scale dependence of the LO cross section with the individual contributions of the partonic 
channels (left panel) and scale dependence of the LO and NLO cross sections (right panel) for 
the pp — > tttt + X process at the LHC for ^/s = 14 TeV. The scale is set to a common value 
M-R = fJ'F = £ • Mo where ^ = 2m ( . 



^-(LHCWeV,™* = 173.2 GeV, MSTW20081o) = 12.056+^!+^ fb (3.1) 

fj^-°(LHCi4TeV,mi = 173.2 GeV, MSTW2008nlo) = fb. (3.2) 

This leaves us with an NLO TC— factor equal to /C = 1.27 and a positive NLO QCD 
correction of the order of 27%. Also given in Table 2 are the integrated NLO cross sections 
for two values of the unphysical cutoff parameter a max that is a common modification of 
subtraction terms in the phase space away from the singularity, first introduced in [58] . To 
be more specific we use two values, a max = 1 that corresponds to the original formulation 
of [56, 57], and 

c^max — 0.01. The independence of the final result on the value of the 
a max parameter is clearly visible in the Table 2. This is a strong consistency check of the 
calculation of the real emission part. For more details on the a max implementation in the 
Helac-Dipoles package see e.g. [22, 43]. 

The theoretical uncertainty of the total cross section, associated with neglected higher 
order terms in the perturbative expansion, can be estimated by varying the renormalization 
and factorization scales in a s and PDFs, up and down by a factor 2 around the central scale 
of the process, i.e. pLQ. An observed change in the value of crr^Q for this two scale choices 
i.e. 0.5^0 an d 2/io is truly asymmetric. Taken very conservatively, as a maximum of these 
two results, the scale uncertainty at LO is estimated to be at the level of 78%. However, 
in case like this it is more appropriate to symmetrize the errors. After symmetrization 
the scale uncertainty at LO is assessed to be instead of the order of 59%. After inclusion 
of the NLO QCD corrections, the scale uncertainty is reduced down to 26%. In Figure 4 
a graphical presentation of the scale dependence is given, both at the LO and NLO. We 
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observe a dramatic reduction of the scale uncertainty while going from LO to NLO. Figure 
4 also shows the scale dependence of the LO cross section with the individual contributions 
of the two partonic channels. At the central scale, /xo = 2mt, the gg channel dominates 
the total pp cross section by about 87% followed by the qq channel with about 13%. 

The theoretical uncertainty as obtained from the scale dependence of the cross section 
is not the only source of systematic uncertainties. Another source of uncertainties comes 
from the parameterization of PDFs. We estimate PDF uncertainties via the so called 
Hessian method where, besides the best fit PDFs, a set of 40 PDF parameterizations is 
provided that describes ±la variation of all parameters which have been used to obtain 
the global fit. These uncertainties are due to experimental errors in the various data that 
are used in the fits. We adopt the prescription from Ref. [61, 62] and based on these 
40 sets calculate asymmetric uncertainties that estimate maximal variations of the result 
calculated for the central value in the positive and negative directions. They amount to 
+5.7% and —4.5%. However, this method does not account for the theoretical assumptions 
that enter into parameterization of the PDFs which are difficult to quantify within a given 
scheme. Therefore, in the next step we use a different PDF set, namely CTEQ PDF [63, 64] 
that should have different theoretical assumptions. Specifically, we employ CT09MC1 
PDFs with 1-loop running a s at LO and CT10 PDFs with 2-loop running a s at NLO. We 
compare the results for the central value using the best fit PDFs. With [iq = 2mt our 
findings can be summarized as follows: 

CT^-°-(LHCi4TeV, m t = 173.2 GeV, CT09MC1) = 11.414(8) fb , (3.3) 

CT^°(LHCi4TeV,m t = 173.2 GeV, CT10) = 14.37(2) fb . (3.4) 

The MSTW2008 results are larger than the CTEQ predictions by 5.6% at LO and 6.7% at 
NLO, which is comparable to the individual estimates of MSTW2008 PDF systematics 2 . 
In fact, since different values of a s are associated with the MSTW2008 and the CTEQ 
sets one should also consider the uncertainties coming from the determination of the value 
of the strong coupling constant, which is fitted together with the PDFs. Clearly, further 
studies are needed to clarify these issues. However, since the PDF uncertainties for the 
process under scrutiny are well below the theoretical uncertainties due to scale dependence, 
which remain the dominant source of the theoretical systematics, we leave such studies for 
the future. 

3.2 Differential cross sections for fiQ = 2mt 

As already mentioned, total cross sections are mostly influenced by final-state production 
relatively close to the threshold as defined by particle masses. On the other hand, differ- 
ential cross sections extend themselves up to energy scales that are much larger than the 
threshold, and may show larger shape distortions in such high-energy regions. Therefore, 
in the next step we turn our attention to the differential cross sections. 

2 Let us point out here that for the old CTEQ6 PDF sets [65, 66] that are still widely used in 
the phenomenological studies at the LHC we have obtained the following results ctJj'Jj(LHCi4TcV, rat = 
173.2 GeV, CTEQ6L1) = 8.259(4) fb, erjgf (LHCi 4Te v , nn = 173.2 GeV, CTEQ6M) = 14.74(2) fb. 



-7- 




800 1200 1600 2000 2400 



H T [GeV] 



Figure 5. Averaged differential cross section distributions as a function of transverse momentum 
of the tt pair (upper-left panel) and the top quark (upper-right panel) for pp — > tttt + X production 
at the LHC with t/s — 14 TcV. Also shown is the differential cross section distribution as a function 
of the total transverse energy of the system (lower panel) . The dash-dotted (blue) curve corresponds 
to the LO, whereas the solid (red) one to the NLO result. The scale choice is fiF = fJ-R = Mo — 2mj. 
The uncertainty bands depict scale variation. The lower panels display the differential K, factor. 

We have checked as many as 16 observables. Here we present three cases where the 
differential JC factor, defined as the bin-by-bin ratio of the NLO result to the LO one for the 
central scale value /j®, has been found to be mostly distorted. In Figure 5 we present the 
averaged transverse momentum distribution of tt pair and of the top quark together with 
the distribution of the total transverse energy of the titi system. The dash-dotted (blue) 
curve corresponds to the LO, whereas the solid (red) one to the NLO result. The upper 
panels show the distributions themselves and additionally include the scale-dependence 
bands obtained with a variation of the central scale by a factor of two. The lower panels 
display the differential /C factor. 
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We observe that at the LHC, employing a fixed scale /xo = 2mj, the NLO corrections 
to the transverse momentum distributions do not simply rescale the LO shapes, but also 
induce distortions at the level of 60%. In case of the total transverse energy of the system, 
which we define as a sum of (anti-)top quark transverse energies 

i=l,2 i=l,2 



E T (t) = yfm 2 t +p T (t), (3.6) 

the situation is even more severe with an observed distortion at the level of 80% or more. 
Clearly, large and negative NLO corrections affect the high Ht region. One can also note 
that the NLO error bands, estimated through scale variation, do not fit nicely within the 
LO ones at lowpr's, as one should expect from a well-behaved perturbative expansion. We 
can summarize our conclusions at this stage by remarking that the fixed-scale choice = 
A*-F = A*o = 2mt does not ensure stable shapes when going from LO to NLO. Therefore, 
differential cross sections are properly described only when the NLO QCD corrections are 
taken into account. 

3.3 Integrated cross section and its scale dependence for fiQ = Ht/^ 

With the goal of stabilizing shapes in the high px and Ht regions, that are relevant for the 
new physics searches, we have explored a dynamical choice for and \xp that helps to 
achieve flatter differential /C-factors, thus describing more appropriately the kinematics of 
the process far away from the threshold. To this end, we have incorporated the dynamical 
scale option into the Helac-Nlo framework. The implementation has been carefully cross- 
checked by reproducing partial and total cross section results for the NLO QCD corrections 
to the pp —> ttbb + X process with a dynamical scale for the setup number I as given in 
Ref. [21]. A per-mille level agreement has been found in all cases. 

For the process at hand, we explored several possibilities and decided in the end to 
consider the dynamical scale hr = fj,p = /jLq = Ht/^- While preserving moderate NLO 
QCD corrections and the dramatic reduction of the theoretical uncertainty, this new scale 
choice turns out to be particularly effective in improving the stability of distribution shapes. 
Table 3 shows the integrated cross sections at LO and NLO for pp — > titi+ X production 
at the LHC with y/s = 14 TeV, using the same settings as before with the only exception of 
the scale choice, this time set to be fio = Ht/4. Our results can be summarized as follows: 

a t L f °-(LHCi4TeV, mt = 173.2 GeV, MSTW20081o) = 13.891 + ^;°^[+^Jj fb , (3.7) 

CT^°(LHCi4TeV, m t = 173.2 GeV, MSTW2008nlo) = 16.87^5Jj^j fb . (3.8) 

The new results are a bit higher, i.e. by 15% at LO and by 10% at NLO compared 
with Table 2, which is perfectly within theoretical error estimates at the corresponding 



3 The sum of the transverse energies of massless outgoing partons and leptons, has already been advocated 
as a good scale choice in the study of NLO QCD corrections to the differential distributions for pp — > V + 3j 



process, where V = W ± ,Z/-y* [34, 35]. 



-9- 




Figure 6. Scale dependence of the LO cross section with the individual contributions of the partonic 
channels (left panel) and scale dependence of the LO and NLO cross sections (right panel) for 
the pp — > tttt + X process at the LHC for ^/s = 14 TeV. The scale is set to a common value 
= = £ • Mo where fx = H T /A. 



Process 


(tlo [fb] 


<Eo =1 [fb] 


^NEo =a01 [fb] 


K-Factor 


[%] 


pp ->■ titt + X 


13.891(9) 


16.87(2) 


16.86(3) 


1.21 


21 



Table 3. Integrated cross section at LO and NLO for pp — > tttt + X production at the LHC with 
y/s = 14 TeV. Results for the MSTW2008 PDF set are presented. In the last two columns the K 
factor, defined as the ratio of the NLO cross section to the respective LO result, and NLO corrections 
in % are given. The scale choice is pip = piR = /.to = Ht/4,. 

perturbative order level. Moreover, the JC— factor obtained with this new scale is smaller, 
of the order of JC = 1.21. The independence of the final result on the value of the a max 
parameter has also been checked in this case as shown in Table 3. For completeness, we 
show in Figure 6 the scale dependence of the LO and NLO cross sections. Again we observe 
a striking reduction of the scale uncertainty while going from LO to NLO. Varying the scale 
up and down by a factor 2 changes the cross section by +80% and —41% in the LO case, 
whereas in the NLO case we obtain a variation of +24% and —25%. With the evaluation 
of NLO QCD corrections the theoretical error has been decreased from about 80% (60%) 
down to 25%. 

3.4 Differential cross sections for fiQ = Ht/^ 

As desired, the results for the integrated cross sections have only slightly changed in com- 
parison with the fixed-scale case. In the following, we study the impact of the different scale 
choice on the differential cross sections. In Figure 7 the averaged differential distribution 
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Figure 7. Averaged differential cross section distributions as a function of transverse momentum 
of the ti pair (upper-left panel) and the top quark (upper-right panel) forpp — > tttt + X production at 
the LHC with y/s = 14 TcV. Also shown is the differential cross section distribution as a function of 
the total transverse energy of the system (lower panel) . The dash-dotted (blue) curve corresponds to 
the LO, whereas the solid (red) one to the NLO result. The scale choice is /j,f — fJ-R = t l o = Ht/^- 
The uncertainty bands depict scale variation. The lower panels display the differential K, factor. 

as a function of the transverse momentum of the ti pair and the top quark together with 
total transverse energy are shown once again, however, this time with the renormalization 
and factorization scales calculated on an event-by-event basis. Instead of 60% distortions 
that we have obtained in the previous case for the pr ti and pr t distributions, the moderate 
and positive corrections of the order or 20% have been established over the whole range 
of pt- An improvement in the differential K— factor is clearly visible in both cases. The 
same conclusions can be drawn for the Ht distribution. In this case the improvement is 
even more impressive. In addition, one can observe that the NLO error bands as calculated 
through scale variation nicely fit within the LO error bands. 
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Figure 8. Averaged differential cross section distributions as a function of transverse momentum 
of the ti pair at LO (left panel) and at NLO (right panel) for pp — > tttt + X production at the LHC 
with yfs = 14 TcV. The dash- dotted (orange) curve corresponds to /io = 2m t , whereas the dashed 
(brown) one to (1q = Ht/^- The lower panels display the ratio of the result with the dynamic scale 
versus the fixed scale. 
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Figure 9. Averaged differential cross section distributions as a function of transverse momentum 
of the top quark at LO (left panel) and at NLO (right panel) for pp — > tttt + X production at the 
LHC with t/s = 14 TcV. The dash-dotted (orange) curve corresponds to /iq = 2m t , whereas the 
dashed (brown) one to /io = -ffr/4. The lower panels display the ratio of the result with the dynamic 
scale versus the fixed scale. 



In general, the idea behind a dynamical scale is to accommodate for multiscale kine- 
matics. In fact, we are not trying to reduce the /C-factor as defined by the ratio of total 
cross sections, but rather to obtain a constant one for the distributions. The fixed order 
approximation is meaningful, when the improved scale choice affects NLO cross sections 
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Figure 10. Differential cross section distributions as a function of total transverse energy of the 
titi system at LO (left panel) and at NLO (right panel) for pp — > titi + X production at the LHC 
with yfs = 14 TcV. The dash- dotted (orange) curve corresponds to fiQ = 2m t , whereas the dashed 
(brown) one to (1q = Ht/^- The lower panels display the ratio of the result with the dynamic scale 
versus the fixed scale. 
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Figure 11. Differential cross section distributions as a function of the invariant mass of the titi 
system (left panel) and the ti pair (right panel) for pp — > titi + X production at the LHC with 
^fs = 14 TcV. The dash-dotted (blue) curve corresponds to the LO, whereas the solid (red) one 
to the NLO result. The scale choice is [If = Mfi = -ffr/4. The uncertainty bands depict scale 
variation. The lower panels display the differential K. factor. 

to a much lower extent than the LO ones. We also demonstrate that with a proper scale, 
we can already obtain good results with a constant rescaling of LO distributions. There 
always remains the question of how these ideas work in practice. A convincing example 
is illustrated on Figure 8 - Figure 10, for three different observables, where fixed scale 
ZC-factors were most unsatisfactory. With the luxury of an NLO cross section, the im- 
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Figure 12. Averaged differential cross section distributions as a function of rapidity of the ti 
pair (left panel) and the top quark (right panel) for pp — > tttt + X production at the LHC with 
y/s = 14 TeV. The dash-dotted (blue) curve corresponds to the LO, whereas the solid (red) one to 
the NLO result. The scale choice is pip = /i/? = fio = Ht/^- The uncertainty bands depict scale 
variation. The lower panels display the differential K, factor. 



provement obtained with a dynamic scale, which can only be guessed at for the lack of a 
NNLO result, is moderate. The change in the shape of LO distributions is, on the other 
hand, rather strong. 

We have compared these two scale choices for many other observables and in all cases 
the new, dynamic one has decreased the shape difference of distributions while going from 
LO to NLO. As an example, the invariant mass distribution of the titi systems and the 
averaged invariant mass distribution of the ti pair are illustrated in Figure 11. Even though 
the invariant mass distributions peak strongly at the corresponding thresholds, there are 
non- negligible tails that are extended to very high invariant mass values. Also for these 
kinematical regions an almost flat differential K— factor has been obtained. In addition, 
in Figure 12, the angular differential distributions are presented. We show in particular 
the averaged rapidity distribution of the ti pair and the averaged rapidity distribution 
of the top quark. As one can see from Figure 12, the ti pairs and the top quarks are 
predominantly produced in the central region. Also here, the differential /C-factor is near 
constant within the whole range of y. 

Through the implementation of the dynamical scale large discrepancies between shapes 
of distributions at NLO and LO have disappeared. New differential K— factors suggest that 
the proper scale choice in LO calculations, which describes the kinematics of the whole 
process on an event-by-even basis, together with a suitably chosen global /C— factor would 
be sufficient for this process. This is good news taking into account that the LO calculations 
are not only less costly, but can be easily merged via CKKW/MLM procedures [67] with 
parton shower programs to obtain complete inclusive hadron level events samples that can 
be directly compared with the experimental data. 
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Figure 13. Differential cross section distributions as a function of the transverse momentum of 
the 1st, 2nd, 3rd and the 4th hardest top quark at LO and NLO for pp — > tttt + X production at 
the LHC with ^fs = 14 TcV. The scale choice is /ip = hr = Ht/^- Also shown is the averaged 
transverse momentum of the top quark at NLO. 



Once the impact of the NLO QCD corrections to the differential cross sections has 
been established, we turn our attention to the properties of the top quarks. In Figure 
13 a comparison of the transverse momentum spectra of the first, second, third and the 
forth hardest top quark in the pp — > titi + X production at LO and NLO is displayed. 
In each plot, the averaged transverse momentum of the top quark at NLO is also shown 
as a benchmark. As expected, the shape of the distributions and their peak change when 
moving from the hardest to the softest top quark configurations. 

Finally, the invariant mass of the two top quarks with the highest pt is given in Figure 
14. This observable is particular interesting in view of new physics searches where the 
mass of a new heavy resonance that decays into the top-anti-top pair is reconstructed as 
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Figure 14. Differential cross section distributions as a function of the invariant mass of the ti 
pair at LO and NLO for y/s = 14 TeV. Left panel: the invariant mass of the two top quarks 
with the highest pt- Right panel: the averaged invariant mass of the ti pair. The scale choice is 
Mf = M-R = H T /4. 



the invariant mass of the two objects with the highest pr in the event. For comparison, 
the averaged invariant mass of the ti pair is also presented. 



4 Summary and Conclusions 

In this paper we have presented a computation of the NLO QCD corrections to four top 
quark production at the LHC. The total cross section and its scale dependence have been 
evaluated for two different scale choices, i.e. for the fixed scale [Ir = [ip = fJ-o = 2> m t and 
for the dynamical scale ixr = = l^o — The impact of the NLO QCD corrections 

on the integrated cross sections is moderate, of the order of 27% for /j,q = 2mt and 21% for 
fj,Q = Ht/4. As to the theoretical uncertainty of our calculation, the contribution related 
to unknown higher-order corrections, as obtained by studying the scale dependence of our 
NLO predictions, is of the order of 25%. We have also analyzed the theoretical error arising 
from different parametrizations of PDFs, being able to quantify it at the level of 5% — 6%, 
thus well below the uncertainty associated with scale dependence. 

Looking only at the total cross section, which is mostly influenced by final state pro- 
duction relatively close to the threshold, both scale choices are in equally good shape and 
the results agree well within the corresponding theoretical errors. On the other hand, dif- 
ferential cross sections show large differences in shape, with distortions up to 80% observed 
within our fixed-scale setting. In particular, large negative corrections are clearly visible 
in the tails of several distributions. Thus, an accurate description of the shapes of observ- 
ables can be given only via full NLO QCD computation in this case. Instead, adopting our 
dynamical scale choice, results have moderate, positive and almost constant corrections of 
the order of 20% for all the investigated observables. This fact suggests that the proposed 
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dynamical scale efficiently accommodates for the multiscale kinematics of the process. 

Well-behaved as it is, the proposed dynamical scale has a number of advantages for 
phenomenological studies. Indeed, it can be used within a LO calculation, together with a 
suitably chosen global K,— factor, to obtain results that well approximate the full NLO QCD 
calculation and can be merged with parton shower programs to obtain realistic hadronic 
events, directly comparable with the experimental data for new physics searches. Partic- 
ularly interesting observables in this sense are the invariant mass and the total transverse 
energy of the ttti system, together with the invariant mass of the two hardest top quarks. 

As a final remark we point out that, despite its relatively small cross section, a good 
theoretical control over the SM pp — > titi background can be phenomenologically relevant. 
Our NLO QCD predictions are of the order of 17 ± 4 [scales] ± 1 [PDF] fb for y/s = 14 
TeV. For comparison, typical predictions of new physics scenarios such as effective four- 
top interactions, Kaluza-Klein gluons or the so-called top-philic Z' (i.e., a Z' which couples 
to third-family quarks only) are set in the range 5 — 100 fb for m new = 1 TeV and 1 — 20 fb 
for m new = 1.5 TeV, where m new is the mass of the new heavy particle or, more generally, 
the energy scale associated with new physics effects. For masses greater than 2 — 3 TeV 
rates are below 1 fb [8]. We believe that a NLO analysis of pp — > titi at the LHC is a 
necessary step towards a correct interpretation of the possible signals of new physics that 
may arise in this channel. 
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